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BOUNDARY  LAYER  FLCW  OF  A  NEWTONIAN  FLUID  SUSPENSION 

Introduction 

In  icebreaking  one  can  consider  three  different  modes^  ramming^  continuous 
or  moving  through  ice  filled  water.  Several  different  continuous  modes  can  be 
identified  depending  upon  the  thickness  of  the  ice.  Analytical  solutions  have 
been  attempted  and  results  obtained  for  the  ships  resistance  in  the  icebreaking 
modes,  however  the  role  of  friction  in  the  resistance  of  a  ship  moving  in  ice  filled 
water  is  not  clear.  It  has  for  the  most  part  been  ignored.  From  practical  exper¬ 
ience  it  is  known  that  high  thrust  is  needed  to  move  a  ship  through  finely  broken 
ice  at  l(3w  speed.  The  reason  for  this  is  not  clear  since  breaking  forces,  buoy¬ 
ant  forces  and  momentum  interchange  would  be  small.  It  was  proposed  that  the 
reason  for  the  high  resistance  was  that  the  finely  divided  ice  can  pack  closely 
together  and  in  essence  change  the  effective  viscosity  of  the  water.  This  is  a 
well  known  effect  of  solid  fluid  mixtures.  Accordingly  some  analytical  efforts  were 

instigated  to  understand  ships  resistance  in  ice  fields  of  various  density  and 
ice  sizes. 

Analytical  Model 

Studies  of  particulate  flows  have  shown  that  when  a  large  number  of  particles 
are  present  in  a  suspension,  a  particle-free  layer  soon  develops  near  the  wall 
(Ref.l).  The  thickness  of  this  region  depends  on  the  mean  particle  concentration  c 
(Fig.III-1),  As  a  result  the  viscosity  of  the  flow  field  close  to  the  wall  is  lower 
than  that  of  the  suspension  far  from  the  wall.  The  same  behavior  has  been  observed 
for  nonspherical  and  spherical  particles  alike  (Ref, 2). 

Figure  Ill-la  indicates  that  the  time  required  to  establish  the  equilibrium 
thickness  is  of  the  order  of  one  minute  or  less  for  high  concentrations. 

Figure  Ill-lb  shows  the  relationship  between  the  equilibrium  layer  thickness 
and  the  concentration.  The  ratio  of  the  particle-free  layer  thickness  to  the 
thickness  of  the  viscous  layer  can  be  approximated  by  the  expression 

<R„) 

The  apparent  viscosity  of  suspensions  is  also  a  function  of  the  particulate 
volume  concentration.  Figure  III-2  (Ref. 3)  shows  the  effect  of  particulate 
concentration  on  the  ratio  of  the  suspension's  apparent  viscosity  to  the 


2. 


C*0  05 


initial 


Final 


C«0  IS 


Initioi 


Find 


Eigure  III-l  The  Development  of  a  Particle-Free  Layer  (of  Thickness  6) 
in  Suspensions  of  Rigid  Spheres.  R  =  0.2  cm^  a  =  0.015  cm 
(a/R  =  0.075),  Q  =  0,356  cm®  sec"^  . 

(a)  Variation  of  6/k  with  time.  The  points  are:  open  circles 
c  =  0.02;  closed  circles:  c  =  0.05;  open  triangles: 

c  =  0.10;  and  closed  triangles:  c  =  0.15, 

(b)  Effect  of  concentration  c  on  60,/^  where  6c»  is 

the  equilibrium  thickness  of  a  particle-free  layer 
attained  after  a  long  time. 

(c)  Tracings  from  photographs  of  the  particle-free  layer  for 
two  suspensions.  After  Karnis,  Goldsmith  &  Mason  (1966b) 
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Figure  III-2  Comparison  of  Theoretical  Viscosity  Relationships  with 
Data  for  Uniform  Spheres 

viscosity  of  the  liquid  phase  alone.  Concentrations  c  ~  0.5  are  considered  to 
be  about  as  large  as  particle  packing  will  permit  before  the  particles  come  into 
intimate  surface  contact  with  each  other.  At  this  limit  a  tenfold  increase  in 
viscosity  is  observed.  The  experimental  values  given  in  Figure  TII-2  are  accu¬ 
rately  correlated  by  the  empirical  expression 

Mw  =  -580  _  (2) 

Equations  of  Motion  for  the  Multilayered  Boundary  Layer 

Figure  III-3  illustrates  the  problem  to  be  analyzed.  The  layer  next  to  the 
boundar3'^  of  thickness  6  is  ice-free.  The  velocity  of  the  interface  between 
the  mush  ice  and  the  ice-free  water  is  r«  The  integral  momentum  equation  for 
this  layer  is  ^ 

I  i  p^6 
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■^6  ' 
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where  is  the  shear  stress  at  the  wall  and  r.,  the  interfacial  shear  stress 
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Figure  III-3 

The  integral  momentuni  equation  for  the  mush-ice  filled  portion  of  the 
boundary  layer  is 

Oy 

P  J  K)  dy'  +  pUv,^  -prv,g  =-T5  (4) 

O 

where  is  the  normal  velocity  component  at  the  outer  edge  of  the  boundary 

layer  and  ,  the  normal  velocity  component  at  the  interface. 

6 

The  normal  velocity  components  at  the  edges  of  the  boundary  layers  can  be 
obtained  from  the  continuity  equations  written  for  each  portion  of  the  boundary 
layer.  The  continuity  equation  for  the  liquid  layer  is 


6 


The  continuity  equation  for  the  ice-filled  portion  of  the  boundary  layer  is 


6 


Combining  Eqs.(5)  and  (6)  yields 


5 


(7) 


Velocity  Profiles 

The  solution  of  the  preceding  integral  equations  requires  the  assumption 
of  plausible  velocity  profiles.  These  profiles  must  satisfy  the  imposed  veloc 
ity  conditions  at  their  extremities  and  the  continuity  of  shear  stress  at  the 
ice-free  water  -  mush  ice  interface. 

If  a  mush-ice  profile  of  the  form 

-j?-  =  ^0  + 

is  assumed  for  laminar  flow  and  the  following  conditions  are  imposed  upon  it 


U«  U  ^  ^  /  T 

=  r;  at  Tl'  =  1 

r  ^ 


dT]' 

then  the  assumed  profile  becomes 


^  =  1  at  =  0 

r 

^  =0  at  n'  =  1 


(9) 


where 


^  =  1  -  +  a^Tl 


-  1  U  =  1  i 
^  "  r  ■  R 


/2 


(10) 

(11) 


The  quantity  R  is  the  ratio  of  the  interfaced  velocity  f  to  the  free-stream 
velocity  U. 

The  ice-free  liquid  velocity  near  the  wall  will  also  be  assumed  to  have 
a  velocity  profile  of  the  form 


^  =  bo  +  +  bgTp  . 


(12) 
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It  is  required  to  meet  the  no-slip  conditions 


u 


r  =  0  at  I)  =  0 


(13) 
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^  =  1  at  71  =  1 


and  also  the  continuity  of  shear  stress  condition  at  the  interface  which  is 

Su, 


lJ.y 


By' 


By 
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=  ^  (-233) 
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=  ^  (b,  +2b3)  . 


(14) 


These  conditions  lead  to  a  laminar  profile  of  the  form 


?  =  ^^  +  (l-b,)Tp 


(15) 


where 


(16) 


The  following  definitions  have  been  incorporated  into  the  previous  equation 

f-' 


U  ^  1 

r  R 

M-v  =  1 

p  G 


(17) 


The  functions  F  and  G  are  known  functions  of  the  mush-ice  volume  concen- 

g 

tration  and  can  be  related  to  Eqs.(l)  and  (2).  For  c  >  0.1^  —  <  0.1^  so  that 

'  <  0.1,  or  9  <  Therefore  for  c  >  0,1  Eq,(2)  becomes 

0  Ov  ^ 


J_  A  =  i(r^ 

5v  Rq  c 


(18) 


7. 


Comparison  of  Eqs.(17)  and  (2)  shows  that 

G  =  .  (19) 

The  ratio  R  will  be  determined  by  the  equations  of  motion  as  a  result  of 
the  analytical  procedures  to  follow^  and  is  therefore  an  implicit  function  of  c 
through  its  dependence  upon  F  and  G. 


Integral  Momentum-Continuity  Equations 

When  the  integral  forms  of  the  momentum  and  continuity  equations  for  the 
ice-free  liquid  layer  are  combined  [Eqs.(3)  and  (5)]  the  resulting  integral 
equation  is 

o'" 

o 

A  similar  but  more  complicated  equation  is  formed  when  the  mush-ice  momentum 

and  continuity  equations  are  combined  [Eqs.(4)  and  (6)] 

1  1 


o'"  I  Xf  -  »  ■‘'l' o'-®  - '■)  S  J  F '’V  pr("  -  D  g  .  T 

0  0 

In  these  equations  the  shear  stress  terms  are  given  by 


"6 


-  -r  (fci) 


(21) 


(22) 

(23) 


If  it  is  further  assumed  that  R  is  independent  of  x  and  that 

1 

5  =  Ax^ 


(24) 


then  Eqs.(20)  and  (21)  reduce  to^  with  the  aid  of  Eqs.(22)  and  (23) , 


J  -  f) '■’1  ■  -0  ^ 


(25) 
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J  r  Vr 


-  dfi'  +  puci -r)a2  J  ^  d-q-puci -R)A® 


2a3,iF 
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The  integrals  which  appear  in  these  two  equations  can  be  evaluated  using 
the  velocity  profiles  (10)  and  (15).  The  assumption  that  R  is  a  function  of  c 
alone  and  independent  of  x  amounts  to  a  closure  of  the  problem  which  leads  to  a 
similarity  solution  for  the  boundary  layer  thickness^  Eq,(24),  For  the  purpose 
of  estiiriating  the  effects  of  the  mush-ice  viscosity^  the  presumption  that  U 
and  r  bear  a  fixed  relationship  to  each  other  depending  on  the  concentration 
alone  does  not  appear  to  be  overly  restrictive. 

The  various  integrals  which  appear  in  the  latter  equations  have  the  follow¬ 


ing  evaluations 


[  H  Hp|  =  1  -t-  W 
0  r  '  3  6 

0 


r  -  fi  -  dll  =  -2-  +  ^ 

J  r  V  r/  ^  15  ^  TS'  3TT 

0 

J  ^  S  -  t) ‘=(5 


With  these  equations  (25)  and  (26)  reduce  to 


pRUA2[i^  +  15  ■  3^]  + 


“g 


LF  I  3  ^  it;  "F 


I  +11  ^  .M. 
r  3j  G 


Eliminating  A  from  the  two  equations  yields 


fb  +  2  £  a  Ti  i  _  W  +  2-l_  F  ^ 

L^i  ^  G  ®  J^F  V  3  ~TT)  G  ^  3j  G  Li5  ^  15  30_. 
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Since  ag  and  are  related  through  the  expression 


Kt-O 


(33) 


Eq.(32)  can  be  written  in  the  form 

bf  +  K.^h  +  Kq  =  0  (34) 

where 

50  F^G  -  20  FG  -  12  G^  +  30  F^ 

4G®  -  10  F®G  - 

and  (35) 

^  _  4F®  -  40  I®G  +  20  FG  +  8G® 

°  4g2  -  10  F^G  -  F® 

Computations 

The  values  of  the  most  significant  parameters  will  be  computed  for  three 
values  of  c  =  0,1^  0,3  and  0,5.  Table  III-l  gives  the  values  of  F  and  G  in 
this  range  as  computed  from  Eqs,(l)  and  (2). 

TABLE  III-l 


c 

F 

G 

R 

B 

0.1 

0.1000 

0.6325 

0.240 

-3.167 

0.432 

0.3 

0.0333 

0.2531 

0.208 

-3.808 

0.223 

0.5 

0.0200 

0.1013 

0.283 

-2.533 

0.206 

The  computed  values  of  b^  which  yield  positive  values  of  R  are  approximately 
unity  for  0.1  <  c  <  0.5.  This  means  that  the  ice-free  water  velocity  profile 
is  essentially  a  straight  line. 

The  evaluation  of  b^  can  be  facilitated  by  the  accurate  approximations 

K,  «-3-5  ^ 

F 

K  2  +  5  -  . 
o  G 


(36) 
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Since  the  coefficient  1.0^  the  ratio  R  can  be  obtained  from  the 


expression 


(37) 


Given  =  1  -  —  and  b^^  1,0  the  value  of  A  can  be  computed  from  either 

Eqs,(30)  or  (31).  An  accurate  approximation  for  A  is  given  by 


The  calculated  values  of  B  given  in  Table  III-l  show  the  expected  trend 
toward  a  thinner  ice-free  water  layer  accompanying  the  increase  in  mush-ice 
concentration.  Since  R  is  a  function  of  G/F  through  Eq,(37)  its  variation  is 
governed  by  the  nonlinear  nature  of  this  ratio  as  shown  in  Table  III-l.  Appa¬ 
rently  the  interfaced  velocity  p  is  not  strongly  influenced  by  changes  in 
remaining  roughly  1/4  of  the  free-stream  velocity  U  over  the  range  of  concen¬ 
trations  investigated. 

One  might  envision  the  role  of  the  ice-free  liquid  layer  next  to  the  wall 
as  essentially  a  ‘'slip  layer"  in  which  the  less  viscous  ice-free  water 
"lubricates"  the  sliding  of  the  more  viscous  mush-ice  over  the  solid  boundary. 

With  this  in  mind  a  comparison  of  the  wall  and  interfaced  shear  stresses 
is  instructive.  Their  ratio  is  given  by 


(39) 


Not  unexpectedly  this  liquid  layer  is  little  influenced  by  inertia  effects^ 
transmitting  the  wall  shear  stress  almost  undiminished  from  the  wall  to  the 
interface.  Practically  the  entire  variation  in  shear  stress  from  to  zero 


u. 


is  accomplished  within  the  mush-ice  portion  of  the  boundary  layer^  and  it  is 
here  too  that  the  major  inertial  effects  occur. 


Total  PraR  Force 

The  shear  stress  at  the  wall  is  given  by  Eq,(23) 


^  |j,RUb^  _  pUx  R  ^ 

~T~  ■  X  n,  B  ^  • 


(40) 


The  total  drag  per  unit  width  of  a  wall  of  length  L  in  the  flow  direction  is 
given  by 

u 

D  =  J  T„  dx  =  ___  =  __  (41) 


where 


pUL 


2go 


pUL  ^^0 


K  = 


^  (l^)] 


2F 


^  #) 


B  .  G 


15 


(42) 


Figure  I1I-4  shows  a  plot  of  K  as  a  function  of  c  corresponding  to  the 

data  given  in  Table  III-l.  The  value  of  K  is  shown  extrapolated  to  the  correct 

limit  for  the  case  of  completely  ice-free  flow  of  K  =  1.3. 

The  preceding  analysis  was  developed  on  the  basis  that  the  flow  was  laminar. 

For  ice-free  flow  this  is  the  situation  when  the  Reynolds  number  =  <  3.5(10^) 

M- 

If  L  =  3,0  ft  and  U  =  6  fps  =  3.55  kts,  ^  10.0(10^)  which  would  lead  one  to 

H 

conclude  that  the  laminar  boundary  layer  is  quite  shorty  which  is  generally  true 
for  single  phase  flow.  It  has  been  shown,  however,  that  turbulence  in  the  free 
stream  greatly  influences  the  transition  to  turbulent  flow  within  the  boundary 
layer  itself.  In  the  case  of  a  ship  the  free  stream  is  actually  at  rest  while 
the  plate  moves  relative  to  it  so  that  free-stream  turbulence  is  negligible. 

This  coupled  with  the  high  viscosity  of  most  of  the  boundary  layer  and  the  fact 


13. 


that  very  little  water  flows  into  the  ice-free  region  could  result  in  laminar- 
ized  boundary  layer  flow  at  Reynolds  numbers  well  above  the  normally  accepted 
critical  value. 

Results  of  Newtonian  Analysis 

The  effect  of  a  thin  layer  of  ice-free  water  separating  the  mush-ice  from 
the  wall  is  to  increase  the  total  drag  relative  to  that  which  is  present  in  the 
complete  absence  of  floating  ice.  As  expected^  this  additional  drag  increases 
as  the  volume  concentration  of  the  raush-ice  increases  to  its  maximum  possible 
concentration.  A  tenfold  increase  in  the  apparent  viscosity  of  the  mush  ice 
is  accompanied  by  an  approximately  four-fold  increase  in  the  total  drag.  This 
significant  increase  has  obvious  implications  for  bodies  moving  through  ice 
floes . 
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NEWTONIAN  FLUID  SUSPENSION 
NOMENCLATURE 

A  =  coefficient,  ice-free  boundary  layer  thickness 
a  =  radius  of  spheres 

ag  =  coefficient,  mush-ice  velocity  profile 
B  =  coefficient,  ice-free  boundary  layer  thickness 
bj  =  coefficient,  ice-free  velocity  profile 
c  =  volume  concentration  of  mush-ice 
D  =  total  drag  per  unit  width  of  plate 
F  =  5/5^ 

G  = 

K  =  coefficient,  total  drag 
Kq  =  dimensionless  coefficient 
=  dimensionless  coefficient 
Q  =  volumetric  flow  rate 

R  =  r/v 

Rq  -  pipe  radius 
U  =  free  stream  velocity 

u  =  x-component  of  velocity,  ice-free  water 

V  =  y-component  of  velocity,  ice-free  water 

Uy  =  x' -component  of  velocity,  mush-ice 

Vy  =  y “^-component  of  velocity,  mush-ice 

v^  y •'-component  of  mush- ice  velocity  at  6^ 

v„  =  V.  =  y^-component  of  inush-ice  velocity  at  6 
0 

v^  ==  v^  =  y ^-velocity  component  at  6,  ice-free  water 

X  =  coordinate  in  flow  direction,  ice-free  water 
=  coordinate  in  flow  direction,  mush-ice 


15. 


y  =  coordinate  normal  to  flow^  ice-free  water 
=  coordinate  normal  to  flow^  mush-ice 
f  =  interface  velocity 

6  =  thickness  of  ice-free  water  boundary  layer 
6^  =  thickness  of  mush-ice  boundary  layer 
11  =  y/5 

n'  =  y'/K 


|j,  =  viscosity  of  ice-free  water 
=  viscosity  of  mush-ice 
=  viscosity  of  suspension  liquid  phase 
p  =  density 

=  shear  stress  at  wall 
=  shear  stress  at  interface 

0 
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NON-NEWTONIAN  SUSPENSION  BOUNDARY  LAYER  FLOW 


Introduction 

In  the  previous  section  the  boundary  layer  flow  of  a  fluid  suspension  (in 
that  case  water  containing  floating  ice  chunks)  was  as  analyzed  by  treating 
the  mixture  as  a  Newtonian  fluid  with  concentration  dependent  viscosity.  The 
results  obtained  from  this  analysis  are  applicable  to  flows  in  which  the  volume 
concentration  of  broken  ice  is  less  than  the  critical  value  of  about  50%^  and 
the  ice  chunks  are  small  compared  to  the  characteristic  length  of  the  boundary. 
These  conditions  could  be  realized  under  certain  circumstances  along  the  sides 
of  the  ship's  hull  somewhat  aft  of  the  bow  section.  In  the  vicinity  of  the  bow^ 
however,  compacting  of  the  ice  by  the  wedged  shape  of  the  bow  could  be  expected 
to  locally  increase  the  concentration  of  the  ice  chunks  to  values  exceeding  the 
critical  concentration,  even  if  the  undisturbed  free-field  concentrations  were 
initially  below  the  critical  value.  There  is,  of  course,  the  possibility  that 
in  practice  undisturbed  free-field  concentrations  exceeding  the  critical  value 
will  also  be  encountered,  in  which  case  the  compacting  effect  on  the  ice  field 
by  the  ship's  motion  would  be  much  less  "localized". 

The  problem  of  ship  motion  in  broken  ice  fields  for  which  the  Newtonian, 
apparent  viscosity  model  is  inappropriate  is  clearly  of  great  interest  and  worthy 
of  further  study.  The  initial  considerations  of  this  difficult  and  apparently 
almost  neglected  problem  will  be  the  subject  of  the  remainder  of  this  report. 

Laboratory  Model  Test 

Before  attempting  any  analytical  investigation  of  highly  concentrated  ice- 
water  suspensions  it  was  decided  that  it  would  be  helpful  if  the  behavior  of 
such  mixtures  was  modeled  and  observed  in  the  laboratory.  A  simple  demonstra¬ 
tion  apparatus  was  therefore  devised.  The  purpose  of  this  experiment  was  to 
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investigate  some  of  the  properties  of  ice  and  water  and  the  effects  of  moving 
a  boat  through  broken  ice.  The  apparatus  used  consisted  of  a  large  tub  5'  by 
2'  by  5'*  deep-filled  halfway  with  water  and  two  small  boats^  the  shorter  made 
of  plastic^  the  longer  of  wood. 

The  first  experiments  were  done  using  crushed  ice.  In  order  to  get  the 
water  in  the  tub  cooled  down  ice  cubes  were  first  allowed  to  melt  in  the  water 
and  then  the  crushed  ice  was  added.  One  of  the  boats  was  then  be  placed  at 
one  end  of  the  tank  and  moved  slowly  through  the  ice. 

Analysis  of  Model  Test 

Instead  of  the  boat  moving  from  one  end  of  the  tank  to  the  other^  breaking 
through  the  ice,  it  was  found  that  the  ice  would  move  in  around  the  boat^  inter¬ 
locking  to  form  an  ice  sheet_,  following  which  the  entire  ice  field  would  move 
with  the  boat.  Even  ice  originally  behind  the  boat  would  cling  to  other  pieces 

of  ice  and  move  along  so  that  the  boat  was  not  only  pushing  the  ice  field  in 

front  of  it  but  pulling  the  one  in  back  of  it  too.  It  appeared  as  if  all  the 

tiny  particles  of  ice  had  frozen  together.  This  was  not  the  case^  however, 
since  when  the  field  was  poked  with  a  finger  all  the  particles  easily  separated 
showing  that  the  chunks  had  not  frozen  together. 

If,  as  the  boat  was  moving,  a  portion  of  the  ice  field  to  the  side  of  the 
boat  was  broken  off  by  running  a  finger  through  the  field  perpendicular  to  the 
side  of  the  boat,  this  separated  part  of  the  ice  field  would  remain  as  a  sheet 
of  ice  and  move  in  behind  the  boat  once  it  had  passed  by. 

With  the  crushed  ice  there  was  never  any  difficulty  reproducing  these  phe¬ 
nomena.  The  best  results  were  obtained  with  a  water  temperature  between  33°F(0-2°C) 
and  38°F  and  a  surface  density  of  607c  or  707o  ice  to  water.  The  same  effects 
were  obtained  with  a  surface  density  of  ice  to  water  as  low  as  507o.  Some  wall 
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effects  were  encountered  during  the  tests.  If  the  ice  was  moved  away  from  the 
wall  at  the  beginning  of  the  boat*s  run^  it  would  form  the  ice  sheet  and  then 
once  again  be  forced  against  the  wall  by  the  motion  of  the  boat. 

The  pictures  on  the  following  pages  illustrate  the  phenomena  described 
above  for  crushed^  tiny  particle  ice.  The  ice  cubes  in  these  photographs  are 
what  remained  of  the  original  cubes  after  cooling  down  the  water.  The  smaller 
pieces  were  subsequently  added  as  crushed  ice.  The  black  spot  in  some  pictures 
are  on  the  bottom  of  the  tub  and  therefore  show  holes  in  the  ice  field. 

Figures  III-6. III-7  -  In  these  pictures  the  boat  started  moving 

from  the  end  of  the  tank.  It  carried  with  it  the  entire  field  of  ice.  This 
cah  be  seen  by  noticing  the  perpendicular  line  directly  behind  the  boat^  indi¬ 
cating  the  end  of  the  ice  field.  The  blur  at  the  bottom  of  the  picture  is  the 
ice  crushing  against  the  opposite  wall. 

Figure  III-8  -  In  this  picture  the  boat  was  placed  in  the  center  of  the 
field  and  then  started  moving.  It  carried  with  it  the  field  in  front  of  it 
and  the  field  in  back  of  it.  The  arrow  points  to  the  line  made  by  the  moving 
field  showing  clearly  how  the  entire  field  moves. 

Figures  III-9  and  IlI-lO  -  Here  the  boat  started  moving  at  the  end  of  the 
tank  carrying  the  entire  field.  A  portion  of  the  field  was  broken  off  by  run¬ 
ning  a  finger  through  the  field  perpendicular  to  the  side  of  the  boat.  One  can 
see  that  as  the  boat  moves  forward  still  carrying  the  field^  this  separated 
portion  moves  in  behind  the  boat.  The  black  spots  are  on  the  bottom  of  the 
tank.  They  show  where  there  are  holes  in  the  field.  The  field  still  moves 
even  with  the  presence  of  a  large  number  of  holes. 


Figure  III-5  Figure  III- 
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Figure  III-7  Figure  III-8 
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Figure  III-9  Figure  III-IO 
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The  crushed  ice  particle  results  looked  almost  identical  for  the  two  dif¬ 
ferent  shaped  boats.  Unfortunately^  very  few  details  or  specific  properties 
of  the  ice-water  flow  could  be  identified  using  such  small  particles.  When 
intermediate  size  ice  cubes  were  used_j  however^  much  more  could  be  seen  and 
the  results  began  to  make  some  sense.  Photographs  of  actual  icebreakers  in 
broken  ice  show  that  there  are  about  50  -  100  pieces  of  ice  to  the  length  of 
the  boat.  This  makes  the  crushed  ice  experiments  somewhat  more  realistic  than 
the  ones  using  ice  cubes  7-15  pieces  to  the  length  of  the  boat. 

With  these  larger  chunks  of  ice  one  could  easily  see  the  “lock-up”  of  all 
the  ice.  Once  again  they  appeared  to  be  frozen  together  but  when  poked  the 
cubes  would  separate.  Even  without  the  boat  in  the  tub  the  ice  would  seem  to 
lock-up  when  a  wave  in  the  water  was  created. 

One  of  the  first  things  investigated  was  whether  or  not  the  particle  size 
had  an  influence,  i.e.,  is  there  a  particle  size  above  which  the  boat  will  not 
push  the  whole  ice  field?  To  test  for  this  effect,  full  size  ice  cubes  were 
initially  placed  in  the  tub.  As  the  cubes  melted  the  boat  was  passed  through 
the  field  of  smaller  and  smaller  ice  chunks.  It  was  concluded  that  if  the 
water  was  cold  enough  (around  33'^F  -  34°F)  and  if  the  cubes  were  in  the  water 
long  enough  to  begin  melting,  the  particle  size  was  of  no  importance.  When  the 
water  was  cold  and  full  size  cubes  were  added  there  was  no  effect  at  first, 
but  given  a  minute  or  two  to  begin  melting,  the  ice  field  could  easily  be  moved 
by  either  boat.  Perhaps  this  melting  creates  a  new  coefficient  of  friction  or 
some  new  resistance  between  the  ice  and  the  water. 

In  order  to  find  the  conditions  which  determined  whether  the  boat  would 
carry  the  field  or  just  break  through  the  ice,  the  boat  was  watched  carefully 
as  it  began  moving.  It  was  observed  that  under  the  correct  conditions  (cold 
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water  and  no  wall  effect  present)  if  the  boat  submerged  a  cube  in  front  of  it^ 
then  it  would  begin  to  move  the  whole  field.  This  action  was  easily  identified 
because  at  this  point  a  small  extra  force  would  be  needed  to  keep  the  boat 
moving  forward.  If  the  boat  were  moved  away  just  enough  to  allow  the  submerged 
cube  to  come  up^  when  the  boat  was  moved  ahead  again  it  would  begin  breaking 
through  the  ice.  Apparently  there  must  be  a  submerged  cube  present  to  cause  the 
boat  to  move  the  field. 

AgaiUj  as  with  the  crushed  ice^  the  boat  was  able  to  move  the  field  to  its 
side  and  behind.  If,  as  before,  a  portion  of  the  field  was  broken  off  to  the 
side,  it  would  move  in  behind  the  boat  remaining  locked  up. 

Figures  III-ll  and  IlX-12  -  One  can  clearly  see  the  arrangement  of  the  ice 
when  it  is  all  '^locked  up".  The  boats  both  started  at  the  end  of  the  tank  and 
have  carried  the  field  with  them. 

Figures  111-13  and  III-14  -  After  the  boat  has  moved  some  distance  within 
the  ice  field,  portions  of  the  field  were  cut  away  along  a  line  running  through 
the  ice  perpendicular  to  the  side  of  the  boat.  Then  as  the  boat  moved  forward 
these  sections  remained  "locked  up"  and  moved  in  behind  the  boat. 

Figure  III-15  -  Here  one  can  see  that  even  with  the  larger  cubes  the  boat 
can  cany  with  it  parts  of  the  ice  field  that  have  formed  behind  it. 

Figures  111-16  and  111-17  -  These  are  close-ups  of  the  boat  in  the  ice 

field.  The  arrow  in  Figure  III-17  points  to  the  cube  that  has  been  submerged. 

Figure  III-18  -  The  boat  has  been  moved  away  but  the  "locked  up"  ice  field 
can  still  be  seen.  The  finger  points  at  the  cube  that  has  a  submerged  cube 

under  it.  The  arrow  points  at  the  submerged  cube  itself. 
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Figure  III-13  Figure  III-14 
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Figure  III-15 
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Figure  III-18 
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Results  of  Model  Tests 

The  large  cubes  probably  do  not  give  a  totally  realistic  view  of  what  would 
happen  in  a  full-scale  ice  field.  However,  they  do  allow  us  to  view  what  we 
cannot  see  with  the  smaller  crushed  ice  particles.  If  these  results  can  be 
scaled  up  to  a  full  size  icebreaker,  they  would  indicate  that  icebreakers  try 
to  drag  entire  ice  fields  with  them.  If  the  ice  is  in  the  process  of  melting, 
the  boat  must  move  the  whole  ice  field.  The  particle  size,  whether  all  large, 
all  Small,  or  mixed,  does  not  seem  to  matter. 

It  may  be  of  some  significance  to  mention  that  the  longer  boat  seemed  to 
be  able  to  submerge  cubes  easier  than  the  shorter  boat  and  was  therefore  more 
likely  to  move  the  entire  ice  field.  Since  broken  off  portions  of  the  ice  field 
did  move  in  behind  the  boat,  maybe  this  accounts  for  the  icebreakers  getting 
stuck  and  not  being  able  to  back  out. 

In  these  experiments  no  effort  was  made  to  keep  the  boats  from  dragging 
ice  fields.  No  data  has  been  taken  to  measure  the  forces  on  the  boats. 

Although  the  extent  of  the  "locked  up"  ice  field  shown  here  may  be  exag¬ 
gerated  relative  to  full-scale  behavior,  these  tests  do  indicate  that  there  is 
a  tendency  of  floating  particulate  material  to  interact  in  such  a  way  as  to 
cause  large  coalesced  regions  to  move  as  a  unit. 

Analytical  Investigation  of  Drag  in  Highly  Concentrated 
Fields  of  Broken  Ice 

In  the  previous  analysis  the  drag  of  a  fluid  containing  a  low  concentra¬ 
tion  of  suspended,  broken  ice  particles  was  investigated.  In  that  case  it  was 
assumed  that  the  presence  of  particulate  ice  would  cause  the  suspension’s  appa¬ 
rent  viscosity  to  increase  depending  on  the  concentration  of  the  ice  particles. 
It  was  demonstrated  that  under  these  conditions  the  drag  on  a  flat  plate  could 
be  substantially  increased  as  a  result  of  the  ice  particles  in  the  suspension. 


30. 


This  earlier  analysis  was  necessarily  limited  to  suspensions  whose  ice 
particle  concentration  is  sufficiently  small^  that  they  experience  little  or 
no  particle-to-particle  interference.  Such  fluids  do  not  tend  to  become 
"structured"  in  the  sense  that  a  solid  is  since  the  particles  do  not  frequently 
come  into  contact^  but  influence  each  others  motion  only  through  the  action  of 
the  intervening  liquid.  In  this  section  a  model  for  highly  concentrated  ice- 
water  suspensions  which  do  tend  to  become  structured  will  be  proposed  and  then 
used  to  obtain  the  characteristics  of  the  flat  plate  boundary  layer  flow  for 
such  a  fluid. 

Rheological  Model 

The  rheological  properties  of  highly  concentrated  suspensions  of  ice  and 
water  do  not  seem  to  have  been  investigated  either  in  the  laboratory  or  the 
field.  However,  even  a  casual  observation  such  as  described  above  of  the 
response  of  concentrated  suspensions  to  the  motion  of  bodies  through  them  indi¬ 
cates  that  they  cannot  be  treated  as  Newtonian  fluids.  Apparently  the  particle- 
to-particle  interference  caused  by  high  particulate  concentrations  causes 
certain  portions  of  the  ice  field  to  behave  more  like  a  solid  than  a  fluid. 

This  effect  is  undoubtedly  caused  by  the  consolidation  of  the  floating  particu¬ 
late  into  a  loosely  interlocking  matrix  in  which  the  water  phase  plays  a  minor 
role  as  the  lubricant. 

Since  a  direct  experimentally  determined  characterization  of  the  stress- 
strain  relationship  for  broken  ice  fields  does  not  seem  to  be  available  at  this 
time  it  was  decided  to  look  for  an  analogous  flow  field  whose  properties  are 
known  and  then  to  model  the  water-ice  suspension  after  it.  The  most  logical 
analogous  flow  is  that  which  occurs  in  fluidized  beds.  Here  the  particulate 
material  is  suspended  in  an  upward  moving  fluid  stream  rather  than  by  buoyancy 
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forces.  Particle-to-particle  interaction  plays  an  important  role  in  the  per¬ 
formance  of  these  beds  and  so  it  is  a  matter  of  some  concern  to  those  who  use 
them.  Of  the  considerable  literature  on  the  subject  of  fluidized  beds  two 
references  proved  to  be  of  particular  interest  with  regard  to  broken  ice  fields 
(Refs. 4^ 7).  In  both  these  references  it  was  demonstrated  experimentally  that 
highly  concentrated  particle-fluid  suspensions  exhibit  non-Newtonian  behavior  to 
a  degree  depending  on  the  concentration  of  the  particulate  material. 

These  suspensions  are  found  to  require  that  a  minimum  shear  stress  be 
exerted  upon  them  before  they  will  yield  continuously  in  the  manner  of  a  fluid. 

In  Ref.4j  which  is  most  applicable  to  the  problem  at  hand^  it  was  shown  that 
the  stress-strain  relationship  for  a  suspension  undergoing  pure  shear  follows 
very  closely  the  classical  constitutive  equation  for  a  Bingham  plastic.  A 
Bingham  plastic  (Fig.III-19)  is  a  material  which  will  not  yield  until  a  minimum 
shear  stress  is  exceeded,  and  then  once  deformation  does  begin  it  follows  a 
linear  stress-strain  rate  relationship  similar  to  a  Newtonian  fluid.  This 
Bingham  plastic- like  behavior  was  particularly  representative  of  those  suspen¬ 
sions  which  contained  the  larger  size  particles. 

It  could  certainly  be  argued  that  for  shear  rate  as  small  as  1.0  sec"”^ , 
which  are  typical  for  ships,  any  suspension  requiring  a  minimum  shear  stress 
could  be  considered  to  be  a  Bingham  plastic.  The  Bingham  plastic  constitutive 
equation  also  happens  to  be  one  of  the  simpler  non-Newtonian  formulations, 
thereby  lending  itself  to  analytical  techniques. 

Thus,  based  on  observations  of  the  behavior  of  an  analogous  material  and 
the  relative  simplicity  of  its  constitutive  equation,  it  was  assumed  as  a  first 
approximation  that  highly  concentrated  ice-water  suspension  obey  the  constitu¬ 
tive  equations  for  a  Bingham  plastic  in  the  range  of  shear  rates  encountered 
in  ship  motion.  The  actual  confirmation  of  the  Bingham  model  and  the  determination 
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Figure  III-19  Bingham  Plastic,  Stress-Strain  Rate 
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Figure  III-20  Region  of  Plastic  Boundary-Layer  Flow  (Unshaded),  with 
Velocity  Profile,  in  the  Moving  Plate  Problem  with 
\  =  0,  The  vertical  scale  is  much  increased  compared 
with  the  horizontal. 
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of  its  two  material  coefficients  must  necessarily  be  left  for  laboratory  con¬ 
firmation  at  some  later  date.  There  is  strong  evidence  from  Refs, 4  and  7  that 
these  coefficients  will  depend  on  the  size^  distribution  and  concentration  of 
the  particulate  ice.  There  is  little  doubt,  however,  that  these  suspensions 
can  be  characterized  as  Bingham  over  some  range  of  shear  rates  in  the  neighbor¬ 
hood  of  zero  once  the  shear  rate-stress  relationship  is  found  from  rheological 
measurements . 

Constitutive  Equations  -  Two-Dimensional  Bingham  Plastic  Flow 

For  one-dimensional  flows  the  equation  relating  the  shear  stress  to  the 
Strain  rate  (see  Fig.III-19)  can  be  expressed  by  the  equations 


There  are  an  infinite  number  of  ways  which  this  special  case  formulation  can  be 
generalized  to  two  and  three-dimensional  flows.  The  most  common  generalization 
used  is  that  proposed  by  von  Mises.  His  equations  are  summarized  in  Ref, 8  and 
are  as  follows: 


P33  =  3KA 

(44) 

Pi'k  =  ;  1  (p/kPi'k)  ^ 

elastic 

(45) 

Pi'k  =  2Tle4  ;  \  (p/„Pi'fc)  >  V®, 

plastic  . 

(46) 

In  the  first  equation  K  is  the  bulk  modulus  of  the  material.  Since  the  pres¬ 
sure  is  defined  as  p  “  "  ^  “  dilatation  =  the  change  in  volume  per 

unit  volume,  the  first  equation  states  that  the  pressure  is  related  to  the 
change  in  volume  of  a  unit  volume  by  the  bulk  modulus . 
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The  components  of  the  strain  tensor  are  given  by 

2  3Xj^/ 

where  are  the  components  of  the  displacement.  The  dilatation  is  given  by 

• 

The  second  equation  applies  in  the  elastic  region  where  the  material  is 
deformed  but  where  the  deformation  rate  =0.  The  modulus  of  rigidity  is 
given  the  symbol  pi.  If  the  stress  components  are  separated  into  two  parts  so 
that  pj^  are  the  deviatoric  components,  then  the  total  stresses  are  given  by 

Pik  Pik  “  P^Ik 

where  is  the  substitution  tensor. 

The  deviatoric  components  of  the  elastic  stress  tensor  (which  are  related 
to  the  modulus  of  rigidity)  are  related  to  the  components  of  the  elastic  strain 
tensor  by  the  second  equation.  This  equation  is  the  standard  constitutive 
equation  of  linear  elasticity  in  which  the  normal  components  contain  the  ef¬ 
fects  of  the  hydrostatic  pressure  p.  This  equation  can  also  be  written  in 
the  form 


Equation  (45)  applies  only  if  a  minimum  state  of  stress  v  is  not  exceeded 
so  that 

i  p/kPi'k  ^  • 

This  limiting  condition  will  be  discussed  more  thoroughly  in  a  succeeding 
paragraph. 

The  third  equation  expresses  the  fact  that  the  material  will  yield  in  a 
viscous  fashion  once  the  critical  stress  condition  v  has  been  exceeded.  The 
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flow  in  this  regime  is  analogous  to  the  flow  of  a  viscous  fluid  where  1]  is  the 
viscosity  (variable  in  this  case)  and 


e  =  i 


The  quantities  are  the  velocity  components  of  the  flow  field. 

One  more  equation  is  needed  in  addition  to  the  first  three^  and  that  is 
the  equation  for  the  variable  viscosity^  comparable  to  that  which  appears  in 
Eq.(42).  The  generalization  proposed  by  von  Mises  is 


T1  =  Til  +  .  /  ,  •  (^7) 

This  form  of  the  variable  viscosity  reduces  to  the  form  of  Eq,(42)  for  one¬ 
dimensional  flow  and  contains  two  material  coefficients^  the  yield  value  and 
the  reciprocal  mobility  7]^ ,  These  are  the  two  quantities  mentioned  earlier 
which  must  be  obtained  experimentally. 


An  equivalent  form  of  the  apparent  viscosity  equation  (47)  is 


Here  it  can  be  seen  that  as  p/kp/k  apparent  viscosity  approaches 

infinity  as  indicated  in  Fig,III-19  for  one-dimensional  flows.  It  can  be  easily 
shown  that  p/kPik  invariant  under  coordinate  rotation  so  that  the  apparent 
viscosity  is  a  proper  scalar.  On  the  other  hand  can  vary  from  point  to 

point  in  the  material  so  that  similar  to  the  one-dimensional  case  the  apparent 
viscosity  in  Eq. (4)  also  varies  from  point  to  point.  Since  ^  the 

generalized  three-dimensional  model  is  further  complicated  by  the  fact  that  the 
variable  viscosity  does  not  simplify  as  neatly  as  the  one-dimensional  version. 
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Boundary  Layer  Solution  for  Slow  Motion  of  a  Bingham  Plastic 

Obviously  the  analytical  solution  of  a  two-dimensional  flow  around  a  body 
shaped  like  the  hull  of  a  ship  is  a  difficult  task  even  when  the  fluid  is  New¬ 
tonian.  The  added  complication  of  a  variable  viscosity,  which  itself  depends 
on  the  kinematics  or  state  of  stress  within  the  flow  field^  probably  precludes 
the  possibility  of  finding  an  analytical  solution  for  flows  around  two-dimensional 
bodies  except  through  the  use  of  a  computer. 

Fortunately  Oldroyd  (Ref. 9)  has  investigated  a  problem  which  has  some 
bearing  on  the  one  at  hand.  The  results  of  his  analysis  will  be  described  here 
since  they  reveal  some  of  the  salient  features  of  the  flow  produced  by  bodies 
moving  through  two-dimensional  fields  of  Bingham  fluids. 

In  this  analysis  the  elastic  regions  will  be  treated  as  essentially  rigid. 

In  this  way  the  elastic  regions  enter  into  the  picture  only  as  boundary  condi¬ 
tions  at  surfaces  of  transition  from  elastic  to  plastic  states  along  what  could 
be  called  a  yield  surface.  In  the  plastic  region  the  flow  will  be  assumed  to 
be  incompressible  so  that 


where  ,  The  components  of  the  strain-rate  tensor  for  the  plastic  region 

are  given  by 


- 

2  / 


For  two-dimensional  flows  the  stream  function  is  defined  by 


u 


^  V  =  ^ 

ay  ’  ax 


e'  =  -  i±-  e'  =  ^ 
**  ax3y  ’  yy  axay 


'  2  V  ay"^ 


so  that 
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Starting  with  equations  of  motion 


5P  + 

ax^ 


Oldroyd  develops  plastic  boundary  layer  equations  of  the  form 


Til 


Sy^ 


+  2J±  ^)\ 

lax  vax  ay/  ay  \ax  ay/f 

^  =  -  2si  A 

ay  ay  Vax  ay/ 


(48) 

(49) 


on  the  assumption  that 

^  «  1.0 
/  yd  \2/3 
\i\vy 

and 

^^1.0 

Til 

where  d  And  V  are  the  characteristic  length  and  velocity  of  the  flow  field. 

Equations  (48)  and  (49)  are  the  equations  for  the  two-dimensional  creeping 
flow  of  a  Bingham  fluid  in  which  all  inertia  terms  are  absent  and  only  the 
largest  stress  related  terms  remain. 

These  equations  were  applied  by  Oldroyd  to  the  solution  of  the  plastic 
boundary  layer  near  the  semi-infinite  thin  plate  shown  in  Fig.III-20.  In  this 
problem  the  thin  plate  moves  with  a  velocity  V  to  the  left  through  an  infinite 
mass  of  a  Bingham  fluid  which  itself  has  zero  velocity  where  undisturbed. 

Equations  (48)  and  (40)  apply  in  the  region  0  <  y  <  6(x),  i.e.^  the  so- 
called  boundary  layer.  The  boundary  conditions  for  the  flow  in  this  region 
are 

u  =  -Vj  v  =  0  on  y  =  0 

u  =  V  =  0,  —  =0  on  y  =  6(x).  (50) 

^  By  j-  V  / 


An  additional  boundary  condition  must  be  imposed  at  the  edge  of  the  boundary 
layer  describing  the  pressure  exerted  by  the  external  elastic  region  on  the 
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edge  of  the  boundary  layer  itself.  This  pressure  will  depend  on  the  stress 
distribution  in  the  elastic  region.  As  a  first  approximation  it  will  be  assumed 
that  the  pressure  is  a  constant,  so  that 

p  =  p3^  (constant)  on  y  =  S(x)  .  (51) 

The  significance  of  this  assumption  with  respect  to  icebreaking  will  be  dis¬ 
cussed  later. 

The  details  of  the  solutions  of  Eqs.(48)  and  (49)  using  conditions  (50) 
and  (51)  are  given  in  Ref. 9  and  will  not  be  repeated  here.  The  results  of  the 
solution  are  as  follows: 


where 


and 


(52) 

u  =  -V(H-2Y)(1  -y)® 

/Tl 

(53) 

(54) 

Y  .  i 

3x  \\vj 

(55) 

/%va/3 

=  Kavx)  • 

(56) 

This  solution  suffers  from  the  same  deficiencies  as  all  boundary  layer  solu¬ 
tions,  As  usual  the  boundary  condition  v  =  0  on  y  =  6(x)  is  not  satisfied 


vx 


exactly.  However,  for  large  values  of 
approximated  since 


v(6)  =  v( 


TliV 

1/3 


this  condition  is  very  closely 


2vx/ 


The  above  solutions  are  also  unreliable  approximations  very  near  the 
plate's  leading  edge  where  the  boundary  layer  approximations  are  always  inac¬ 
curate.  The  region  of  the  flow  field  in  which  the  approximations  fail  are 
fortunately  quite  small,  and  generally  do  not  prove  to  be  troublesome  if  one 
is  only  interested  in  the  global  features  of  the  flow  such  as  the  total  drag 
on  the  whole  plate. 
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A  curious  feature  of  this  flow  field  solution  is  the  fact  that  the  velocity- 
profile  in  the  boundary  layer  yields 


sy 


y^O 


=  0 


so  that  the  shear  stress  on  the  boundary  is  simply  \)  and  the  total  drag  on  the 
plate  (both  sides)  is  2v  times  the  area  of  the  plate  (one  side) .  This  result 
is  a  consequence  of  the  assumption  of  slow  flow.  It  also  suggests  a  means  by 
which  the  material  coefficient  v  might  actually  be  measured  in  the  laboratory. 

Although  the  flow  around  a  thin  plate  does  not  model  the  flow  around  the 
hull  of  a  ship  it  at  least  indicates  the  nature  of  the  effects  which  the  rhe¬ 
ology  of  the  material  has  on  the  drag.  For  this  simple  configuration  we  find 
that  the  drag  is  dependent  upon  a  property  of  the  suspension  v,  the  area  of 
the  surface^  and  independent  of  the  normal  stress  and  the  velocity  of  the 
plate.  This  leads  to  a  friction  law  of  the  type  F  cc  v  X  area  which  applies  once 
motion  occurs.  To  initiate  the  motion  from  rest  one  must^  of  course^  apply  a 
force  equal  to  v  X  area  otherwise  the  entire  ice-water  field  is  deformed  elastic- 
ally  while  the  plate  displaces  only  slightly  as  it  loads  the  field.  The  fric¬ 
tional  behavior  here  is  analogous  to  the  Coulomb  friction  case  where  sliding 
is  initiated  only  if  a  certain  minimum  tangential  force  is  applied. 

The  "sticking  effect"  caused  by  v  has  particular  significance  for  ship 
hulls  since  it  indicates  that  a  certain  minimum  thrust  must  be  exerted  by  the 
propellers  to  set  the  ship  in  motion  when  it  is  at  rest  in  a  concentrated 
broken  ice-water  field.  This  is  contrary  to  the  behavior  expected  in  less  con¬ 
centrated  broken  ice  fields  where  even  the  slightest  thrust  will  set  the  ship 
into  motion  through  the  Newtonian  fluid. 

The  solutions  given  above  suffer  from  one  other  deficiency  and  that  is  the 
fact  that  the  pressure  on  the  plate  is  found  to  be  less  than  the  pressure  in 
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the  free  stream, 
is  given  by 


When  Eq. (52)  is  evaluated  at  y  =  0  the  pressure  on  the  plate 


This  means  that  the  free  surface  of  the  fluid  would  tend  to  sink  down  as  the 
plate  surface  is  approached.  This  is  certainly  contrary  to  the  observed  be¬ 
havior  along  two-dimensional  bodies  where  there  is  a  definite  "ridging”  upward 
of  the  ice^  especially  near  the  leading  edge.  Oldroyd  (Ref. 9)  shows  that  an¬ 
other  pressure  boundary  condition  at  the  edge  of  the  boundary  layer  is  admis¬ 
sible  which  is 

P=Pi  P2(^)  on  y  =  6(x)  (57) 

where 

P2  -  2Xv  ^ 

X  =  a  dimensionless  constant. 


With  this  modification  the  previous  solutions  are  amended  to  read: 


5(x)  =  3x^' 


TlxV 


2(1  +  3X)vx 
T1,V  \l/3 


-  Pi  = 


1/3 


(1  -  X  -  Y) 


V(H-3X  +  2YU1  -Y)^ 
1+3X 


"  ■  TTii:  '  1 

y  ^  y_  /2(1  +  3X)vx 

3x  V  TliV 


(58) 

(59) 

(60) 

(61) 

(62) 


The  condition  that  and  5(x)  are  positive  requires  X  to  be  a  positive  number 

oy 

The  total  drag  on  the  plate  (both  sides)  in  this  case  now  becomes 

'HiV 


F  =  2vx|i  +  12X^- 


2(1  +  3X)VX 


)n . 


(63) 


Under  these  circumstances  the  value  of  pg (x)  at  the  outer  edge  of  the  boundary 
layer  is  given  by 
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Ps  (x) 


/  TliV 

4Xs)^2(i  +  3X).^x 


1/3 


(64) 


This  result  shows  that  the  pressure  at  the  outer  edge  of  the  boundary  layer  is 
quite  large  near  the  leading  edge  of  the  plate,  and  approaches  a  constant  value 
as  X  becomes  large. 

The  pressure  at  the  plate  surface  is  now  given  by 

p(x,0)  =p,  +4^,(X-1)(2(^T3^)  .  (65) 


The  pressure  at  the  plate  is  once  again  less  than  the  pressure  at  the  outer 
edge  of  the  boundary  layer,  which  in  turn  is  higher  than  the  far  field  pres¬ 
sure  pj_ If  X  >  1,  "ridging"  occurs  throughout  the  boundary  layer.  The  high¬ 
est  portion  of  the  ridge  occurs  at  the  outer  edge  of  the  boundary  where  there 
is  a  rather  abrupt  return  to  the  level  of  the  plastic  region  along  the  yield 
boundary.  Considering  the  thinness  of  the  boundary  layer  itself  it  would 
appear  that  the  nonuniforra  stressing  of  the  elastic  portion  of  the  ice-water 
field  would  probably  cause  the  ridging  to  be  confined  to  a  narrow  region  adja¬ 
cent  to  the  plate  as  expected. 

It  is  also  noted  that  in  this  case  that  the  drag  has  a  modest  velocity 
dependence  as  well  as  a  minimum  value  for  incipient  motion.  This  is,  of  course, 
a  somewhat  more  realistic  result  from  an  intuitive  point  of  view. 


Conclusions  and  Recommendations 

1.  A  search  of  the  existing  literature  on  ice-water  suspensions  indicates 
that  the  rheology  of  such  suspensions  has  not  been  investigated.  A  useful 
characterization  of  highly  concentrated  broken  ice  flow  fields  cannot  be  ac¬ 
complished  until  this  is  done. 
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2.  In  the  absence  of  a  rheological  description  of  broken  ice  fields  it  is 
hypothesized  that  these  fields  will  behave  as  a  Bingham  plastic^  at  least  for 
some  range  of  shear  rates  near  zero.  This  hypothesis  is  based  on  the  examina¬ 
tion  of  analogous  flow  fields. 

3.  The  theoretical  basis  of  two-dimensional  Bingham  plastic  flows  is  found 
to  be  well  established  and  could  be  used  to  determine  the  drag  on  a  ship  hull 
in  a  concentrated  broken  ice  field.  The  complexity  of  the  rheological  model 
and  the  flow  field  geometry  indicate  that  the  only  practical  way  to  carry  out 
such  an  investigation  would  be  to  use  a  computer. 

4.  The  analytical  solution  of  the  slow  flow  around  a  flat  plate  for  a 
Bingham  plastic  reveals  a  mechanism  whereby  a  "sticking  effect"  between  the 
ice-water  suspension  and  the  hull  can  be  explained, 

5.  The  flat  plate  boundary  layer  solution  also  yields  a  hydrostatic 
pressure  variation  within  the  boundary  layer  which  can  be  used  to  explain  the 
phenomena  of  "ridging"  near  the  solid  surface. 

6.  The  fact  that  the  quantitative  results  of  the  flat  plate  boundary  layer 
analysis  are  qualitatively  in  reasonable  agreement  with  observations  indicates 
that  the  Bingham  plastic  model  warrants  further  consideration  as  the  rheological 
model  for  ice-water  suspension  flows.  What  remains  to  be  done  is  to  determine 
the  material  coefficients  \  and  v  and  to  apply  the  two-dimensional  equations 

to  more  realistic  geometries. 
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NOMENCLATURE 

d  “  characteristic  length 
e^^  =  component  of  strain  rate  tensor 
K  =  bulk  modulus 
i  =  plate  length 
P^Pi  “  pressure 

=  component  of  the  stress  tensor 
Uj  =  component  of  strain 
u^v  =  velocity  components 

V  =  plate  velocity 
x^y  =  coordinates 

5  =  boundary  layer  thickness 
A  =  dilatation 

component  of  strain  tensor 
T\  =  apparent  viscosity 
-  reciprocal  mobility 
\  =  constant 
ji  =  modulus  of  rigidity 

V  =  yield  value 

\|;  =  stream  function 


